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ABSTRACT 

Ultra-high spectral resolution observations of time-varying interstellar absorption to- 
wards K Vel are reported, using the Ultra-High Resolution Facility on the Anglo- 
Australian Telescope. Detections of interstellar Cai, Can, Ki, Nai and CH are ob- 
tained, whilst an upper limit on the column density is reported for C2. The results 
show continued increases in column densities of K i and Ca i since observations ~ 4 yr 
earlier, as the transverse motion of the star carried it ~ 10 AU perpendicular to the line 
of sight. Line profile models are fitted to the spectra and two main narrow components 
(A & B) are identified for all species except CH. The column density A^(Ki) is found 
to have increased by 82tg°% between 1994 and 2006, whilst Af(Cai) is found to have 
increased by 32 ±5% over the shorter period of 2002-2006. The line widths are used to 
constrain the kinetic temperature to Tk,K < &Tlt)iyV^ and Tk,B < 114^]^4K. Electron 
densities are determined from the Cal/Call ratio, which in turn place lower limits 
on the total number density of tt-a ^ 7 x lO'^cm"'^ and ne > 2 x lO^cm"^. Calcium 
depletions are estimated from the Cai/Ki ratio. Comparison with the chemical mod- 
els of Bell et al. (2005) confirms the high number density, with n = 5 x 10"'cm~'^ for 
the best-fitting model. The first measurements of diffuse interstellar bands (DIBs) to- 
wards this star are made at two epochs, but only an upper limit of < 40% is placed on 
their variation over 9 years. The DIBs are unusually weak for the measured E(^b-v) 
and appear to exhibit similar behaviour to that seen in Orion. The ratio of equivalent 
widths of the A5780 to A5797 DIBs is amongst the highest known, which may indicate 
that the carrier of A5797 is more sensitive to UV radiation than to local density. 

Key words: ISM: individual objects: k Vel cloud - ISM: structure - ISM: atoms - 
ISM: molecules - ISM: lines and bands - local interstellar matter 



1 INTRODUCTION 

Over the past two decades an increasing body of ob- 
servational evidence has indicated that the diffuse inter- 
stellar medium (ISM) is highly structured on scales be- 
low one parsec, contrary to earlier theoretical expectations 
(McKee & Ostriker 1977) and observations which had sug- 
gested an apparent minimum cloud size (e.g. Crovisier et al. 
1985). Evidence for small-scale structure (SSS) has come 

* E-mail: ktsOras. org.uk; present address: Royal Astronomical 
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from observations of interstellar absorption at radio, optical 
and ultraviolet wavelengths along closely separated lines- 
of-sight. Such sight-lines have been provided by spatially- 
resolved observations of extended extragalactic radio 
sources (Dieter et al. 1976; Diamond et al. 1989; Lazio et al. 
2009), binary stars (Meyer 1990; Watson & Meyer 1996; 
Cordiner et al. 2006) and star clusters (Andrews et al. 2001; 
van Loon et al. 2009). Alternatively, a single background 
source can be observed repeatedly over several years, with 
each epoch probing a slightly different line-of-sight due 
to the transverse proper motion of the source and/or the 
intervening interstellar material. This technique has been 
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applied to radio observations towards pulsars (Frail et al. 
1994; Stanimirovic et al. 2003, 2010) and optical observa- 
tions towards nearby early-type stars (e.g. Blades et al. 
1997; Price et al. 2000; Banks et al. 2001; Welty 2007, see 
also the review by Crawford 2003). These observations in- 
dicate that large differences in column density can exist be- 
tween lines-of-sight with transverse separations from as few 
as tens to thousands of AU. 

Heiles (1997) considered possible theoretical scenarios 
to explain these variations without grossly violating pressure 
equilibrium with the surrounding ISM. His preferred model 
consisted of a population of cold and relatively dense struc- 
tures in the form of sheets or filaments; the observed varia- 
tions could then be produced if a structure was present and 
closely aligned to the line-of-sight. The model required very 
low temperature gas (~ 15 K) to maintain pressure equilib- 
rium. However, evidence for such a low temperature has not 
been found; instead more recent observations have indicated 
that a small but significant fraction of diffuse gas is at high 
pressure (Jenkins & Tripp 2001, 2011). Heiles & Stinebring 
(2007) have summarised the observational and theoreti- 
cal constraints on models for SSS, and theoretical atten- 
tion has now shifted towards solutions which do not re- 
quire dynamical equilibrium (e.g. Heiles 2007). Such over- 
pressured clouds are expected to evaporate quickly, so a 
mechanism (e.g. magnetohydrodynamic waves, as discussed 
by Hartquist et al. 2003) must be invoked to replenish the 
population of small-scale structures. 

The existence of small-scale relatively dense regions 
within the diffuse ISM has important implications for inter- 
stellar chemistry. Although the interstellar radiation field is 
expected to be essentially unattenuated within such struc- 
tures, the high density can result in enhancement of the 
molecular content (Bell et al. 2005). This provides a poten- 
tial method for obtaining new information on the phys- 
ical conditions, and a possible solution to some of the 
challenges which exist in reproducing the abundances of 
molecules observed along lines-of-sight through the diff'use 
ISM (Cecchi-Pestellini et al. 2009). A comprehensive under- 
standing of the physical conditions within the SSS is re- 
quired, so that it may be fed into such chemical models. 

This paper is organised as follows: in section 2 we de- 
scribe previous studies of SSS towards hi Vel, whilst our ob- 
servations are detailed in section 3. The results are presented 
in section 4, along with a line profile analysis. Section 5 con- 
tains a comparison with earlier observations and a discus- 
sion of the physical and chemical conditions towards k Vel. 
Conclusions are presented in section 6. 



2 BACKGROUND 

K Velorum (HD 81188) is a bright V = 2.46 single-lined 
spectroscopic binary star with an orbital amplitude of 
46.5 kms"^ and period of 117 days (Pourbaix et al. 2004); 
the spectral type of the primary is B2 IV, whilst the na- 
ture of the secondary is unknown. The distance to the sys- 
tem obtained from its trigonometric parallax measured by 
the Hipparcos satellite is 165 ± 13 pc (Perryman 1997)^; the 

^ van Leeuwen (2007) presented an updated reduction of the Hip- 
parcos data with significantly improved astrometry for most stars. 



reddening is E^g-v) = 0.10 (Cha et al. 2000). Early high- 
resolution observations of the ISM towards this star taken 
in 1989 and 1994 have been presented by Crawford (1991) 
and Dunkin & Crawford (1999) respectively. 

Small-scale structure towards k Vel was first noted by 
Crawford et al. (2000), who serendipitously discovered that 
the equivalent width of the interstellar Kl line had in- 
creased by ~ 40 per cent between 1994 and 2000, with a 
corresponding but smaller increase in absorption by Nal. 
The enhancement was found to be in a single narrow {b = 
0.24 km s~^) absorption component located at a velocity"^ 
vq — -1-8.5 km s~^. The observations were interpreted as 
indicating that as the proper motion of n Vel carried it 
~ 15 AU tangentially between 1994 and 2000 the line-of 
sight gradually intercepted a filament of cold and/or dense 
material aligned along the line-of-sight, of the type proposed 
by Heiles (1997). 

Follow-up observations by Crawford (2002) detected 
Cal and CH absorption at the same velocity, and found a 
continued increeise in the Kl column density. From a con- 
sideration of the Cal/Call ratio, CH abundance and stellar 
proper motion, Crawford (2002) concluded that the absorb- 
ing material lay in a cool (r~ 100 K), dense (nn > 10^ cm"^) 
cloud with approximate dimensions of 100-1000 AU in depth 
(along the line of sight) and 15 AU in the transverse direc- 
tion. 

Consideration of this cloud prompted Bell et al. (2005) 
to develop a chemical model of such a filament. The column 
densities of several atomic and molecular species were calcu- 
lated for a cloud with physical extent similar to that inferred 
by Crawford (2002) as a function of time, radiation field and 
density. The models which most closely reproduced the ob- 
servations towards n Vel were for a young (~ 50 yr), tran- 
sient (lifetime < 100 yr), and high density (nn > lO^cm""^) 
filament immersed in a typical ambient interstellar radiation 
field. The models predicted that detectable column densities 
of additional diatomic molecules could be present, including 
C2 and OH. 

The line of sight towards k Vel offers an excellent op- 
portunity to explore the on-going chemical and dynamical 
evolution of small-scale structure. We therefore set out to 
re-observe the chemical species studied by Crawford (2002), 
and to conduct a search for C2, in order to further con- 
strain the physical and chemical conditions and test the 
chemical model of Bell et al. (2005). In particular, detec- 
tion of the homonuclear molecule C2 could be used to 
measure the local kinetic temperature through the popula- 
tions of each rotational J-state (and radiation temperature 
and density, if sufficient high- J lines were also detected, cf. 
van Dishoeck & Black 1982). The resolved intrinsic profiles 
of the narrow interstellar lines could then be used to infer the 
turbulent contribution to the line widths (Crawford 1997), 



However, in the case of k Vel the new reduction resulted in a 
'stochastic' solution as opposed to the fully converged solution 
found in 1997. The distances and proper motions calculated from 
the two reductions are consistent to within their mutual error 
bars. For this particular star we have therefore chosen to adopt 
the astrometry from Perryman (1997). 

^ All velocities quoted in this paper are heliocentric. To convert 
to LSR velocities towards ft Vel subtract 11. 43 km s"-'^. 
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Table 1. Summary of the UHRF observations towards k Vel. Exp 
time is the total exposure time, S/N is the resulting per-pixcl 
continuum signal-to-noise ratio and W is the equivalent width 
integrated over the entire velocity profile. Errors are la. The C2 
upper limit is Scr on each of the five lines which fall in the observed 
wavelength region. 





UT date 


Exp time 


S/N 


W 






(hr) 




(mA) 


Cal 


2006 March 10 


4.0 


380 


1.13 ±0.02 


Call 


2006 March 11 


0.5 


160 


23.24 ±0.26 


Kl 


2006 March 09,10 


2.75 


210 


7.49 ± 0.34 


Nal 


2006 March 11 


0.67 


190 


60.18 ±0.14 


CH 


2006 March 09,11 


9.0 


600 


0.25 ±0.03 


C2 


2006 March 12-15 


15.0 


800 


< 0.11 



which could be particularly important for understanding the 
origin and evolution of the small-scale structure. 



3 OBSERVATIONS 

Observations of interstellar absorption towards k Vel were 
made on the nights of 2006 March 9-15 using the Ultra- 
High Resolution Facility (UHRF, Diego et al. 1995) located 
at the coude focus of the 3.9 m Anglo- Australian Telescope 
(A AT). The UHRF is an echelle spectrograph which is ca- 
pable of ultra-high resolution but can only record simul- 
taneous data over a very narrow wavelength region, so a 
separate setup is required for each species. A summary of 
the observations is presented in table 1. New observations 
of the absorption lines due to Ca I, Ca 11 (K line) , K I (7698 A 
line) and Nal (Di line) were obtained, along with rotational 
lines of the A^A ^ X^H (0,0) band of the CH molecule in 
the R2(l) A-doublet (comprising the R/e(l/2) and R//(l/2) 
lines). A search was performed for the C2 molecule in sev- 
eral low-J lines (R(0), Q(2), Q(4), Q(6) and P(2)) of the 
A^n„ ^ X^E+ Phillips (2,0) band, which were covered in a 
single instrumental set-up. 

The detector used for the atomic lines and CH was the 
blue-optimised EEV2 CCD (2048 x 4096 13.5-nm pixels). 
The EEV2 was operated in normal readout mode, which 
provides an RMS readout noise of 3.2 e~ /pixel and an un- 
binned readout time of 180 seconds. Because the image sheer 
(see below) projects the spectrum onto ~ 1400 unbinned 
pixels, the EEV2 was binned by a factor of eight perpen- 
dicular to the dispersion to reduce the readout noise and 
time. The EEV2 CCD suffers from severe fringing at wave- 
lengths above ~ 6500 A which particularly affected the K I 
observations, necessitating careful flat-fielding. For the C2 
observations the red-optimised MITLL3 CCD (2048 x 4096 
15- vim pixels) was used due to its higher quantum efficiency 
and lower fringing at long wavelengths. As the MITLL3 de- 
tector could not be binned it was operated in slow readout 
mode, which provides a readout noise of 1.8e~/pixel and a 
readout time of 300 s. 

The UHRF was operated in 'R= 1E6' mode, the highest 
spectral resolution available, which is vital to study the pro- 
files of narrow interstellar lines. Wavelength calibration used 
a Th-Ar arc lamp, whilst flat flelds were obtained with an 
internal quartz-halogen calibration lamp. In this resolution 
mode the instrument resolves the lines from the arc lamp, so 
a stabilised He-Ne laser was used to measure the instrumen- 



tal resolution. The resolving power {R = A/AA) was found 
to be i? ~ 845, 000 for the EEV2 and R ~ 940, 000 for the 
MITLL3, corresponding to velocity resolutions of 0.355 and 
0.319 km s^^ respectively. 

Individual exposures were limited to 30 min to prevent 
degradation of the instrumental resolution by the chang- 
ing heliocentric velocity correction due to the rotation of 
the Earth. The detector readout times therefore constituted 
a significant overhead, and detector readout noise made a 
small but significant contribution to the overall noise level. 

The UHRF incorporates an image sheer which refor- 
mats a 1.5" X 1.5" entrance aperture to a 0.06" slit to re- 
duce slit losses, but itself has a peak throughput of only 
15%. Observing conditions were clear on the nights of 9 and 
15 March, but those of 10-14 March were frequently inter- 
rupted by cloud and fog. The seeing during the observations 
varied in the range 1-3", mostly around 1.5". 

Data reduction was performed using the noao . onedspec 
routines in iraf'^. Pixel-to-pixel sensitivity variations and 
fringing were removed by division by a normalised flat- 
field. The scattered light background was measured from the 
inter-order regions and subtracted. The Nal, Cal, Kl, CH 
and C2 observations were corrected for telluric absorption 
by division by a well-exposed spectrum of a bright, fast- 
rotating, early- type standard star {a Vir and/or a Eri); the 
Ca 11 observations were found to be free of telluric contam- 
ination. Individual exposures were converted to a heliocen- 
tric reference frame before being co-added and converted to 
a velocity scale. 

Additional observations of diffuse interstellar bands 
(DIBs - see the review by Sarre 2006) were obtained us- 
ing other spectrographs at two additional epochs. The ear- 
lier observations were obtained during January 1995 as part 
of a survey of DIBs towards Southern stars conducted by 
one of us (SJF). The instrument used was the coude echelle 
spectrograph on the 74-inch telescope at the Mount Stromlo 
Observatory. The resolving power was 7? ~ 60, 000 and the 
signal-to-noise was ~ 550. Later observations were obtained 
during June 2004, as part of a larger study of DIBs as tracers 
of SSS (Cordiner et al. 2006, Cordiner et al. in prep). The 
instrument used was the UCLES spectrograph on the AAT. 
The spectral resolution was i? ~ 58, 000 and the signal- 
to-noise was ~ 2, 300. Full details of the observations and 
data reduction have been given by Cordiner (2006) and 
Cordiner et al. (2006). The 1995 and 2004 data were cor- 
rected for telluric absorption using composite spectra ob- 
tained by co-adding spectra of the same two stars (a Vir and 
Q Eri), observed at both epochs. These telluric standards 
were obtained at similar airmasses to k Vel and required 
negligible intensity scaling, so any residual stellar features 
should be the same at both epochs. 



4 RESULTS & LINE PROFILE ANALYSIS 

High signal-to- noise (S/N) was achieved in the UHRF obser- 
vations of each line, improving upon the earlier observations 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Interstellar Cal absorption observed towards k Vel. In the upper panel, the thin histogram is the observed spectrum and the 
thick lino is the fitted model profile. The lower panel shows the residuals. The dashed line indicates the continuum level and tick marks 
indicate the positions of each absorption component listed in table 3. The velocity scale is heliocentric. 
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Figure 2. As in figure 1, but for interstellar Call K absorption. 
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Figure 3. As in figure 1, but for interstellar Nal Di absorption. The model does not completely fit the core of the line, see section 4.2.3. 
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Figure 4. As in figure 1, but for interstellar Kl absorption. 
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Figure 5. As in figure 1, but for interstellar CH absorption. 



of Crawford (2002) by factors of 1.5-2 in S/N. All of the 
atomic absorption lines were detected at high significance. 
Absorption by the CH molecule was also found but none 
of the C2 hnes were detected, for which we determined up- 
per limits. The UHRF spectra are shown in figures 1-5; the 
continuum signal-to-noise ratios and measured total equiv- 
alent widths are given in table 1. Preliminary results from 
these observations were presented by Smith (2010); they are 
superseded by the results presented here. 

The line profiles of the detected species were mod- 
elled using the vapid interstellar line modelling program 
(Howarth et al. 2002), which computes the optimum col- 
umn density A'^, velocity v and Doppler broadening param- 
eter* h for a given set of absorption components using a 
non-linear minimisation algorithm. Individual absorption 
components have a Voigt profile formed by the convolution 
of the Lorentzian natural line shape, the instrumental point- 
spread-function (taken to be Gaussian with widths given by 



* The broadening parameter b is related to the FWHM of the 
corresponding Gaussian distribution by fe — F^HM 



2Vlri2 



the velocity resolution found in section 3) and a Gaussian 
with broadening parameter b. 

The wavelengths and oscillator strengths adopted in the 
analysis are given in table 2, and were taken from Morton 
(2003) for the atomic species, Black & van Dishoeck (1988) 
for CH and Sonnentrucker et al. (2007) for C2. The CH ro- 
tational lines form a A-doublet split by 1.43 km s~^; the ve- 
locities output by VAPID are with respect to the bluest tran- 
sition, and were converted to the weighted mean of the two 
lines by subtraction of 0.72 kms~^. The two levels which 
contribute to the doublet were assumed to have equal popu- 
lations, because the energy splitting between them is negligi- 
bly small compared to the thermal energy (cf. Crane et al. 
1995); the oscillator strengths used in vapid were appor- 
tioned accordingly. The Kl and Nal Di transitions are split 
by nuclear hyperfine structure, resulting in two lines 0.370 
and 1.083 km s~^ apart respectively. The relative strengths 
of the transitions are equal to 2F+1, which for both observed 
Nal and Kl transitions gives a ratio of 5:3 with the stronger 
transition being the redder of the pair (see Welty et al. 1994 
or Morton 2003 for derivations). The VAPID output veloci- 
ties were converted to the weighted mean velocity of the 
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Table 2. Atomic and molecular line data adopted in the line 
profile analysis. Wavelengths and oscillator strengths / are taken 
from Morton (2003) for the atomic species, Black & van Dishoeek 
(1988) for CH and Sonnentrucker et al. (2007) for C2. The Nal 
and Kl lines have hyperfino structure, whilst C2 and CH each 
have several rotational lines. Oscillator strengths have been ap- 
portioned between hyperfine components in a 5:3 ratio according 
to the prescriptions of Welty et al. (1994). For C2, the line ratios 
given by Sonnentrucker et al. (2007) have been adopted. 







Wavelength 

(A) 


/ 


Cal 




4226.728 


1.77 


Call 


K 


3933.6614 


6.267 X 10-1 


Ki 




7698.9586 
7698.9681 


1.248 X 10^1 
2.079 X 10^1 


Nal 


Di 


5895.9109 
5895.9322 


1.200 X 10-1 
2.001 X 10-1 


CH 


%/(V2) 

R/e(l/2) 


4300.3235 
4300.3030 


2.53 X 10-3 
2.53 X 10-3 


C2 


R(0) 
Q(2) 
Q(4) 
Q(6) 
P(2) 


8757.686 
8761.194 
8763.751 
8767.759 
8766.031 


1.40 X 10-3 
7.0 X 10-* 
7.0 X 10-* 
7.0 X 10-* 
1.4 X 10-* 



transitions by subtraction of 0.231 kms- and 0.677 km s- 
for Kl and Nal respectively. For each C2 J-state, the os- 
cillator strength was apportioned between the P, Q and R 
branches according to the relative line strengths given by 
Sonnentrucker et al. (2007). The results of the fitting are 
presented in table 3. 

4.1 Uncertainties 

The results of the line profile analysis are affected by a num- 
ber of sources of uncertainty, both statistical and system- 
atic. Rigorous estimates of the statistical uncertainty were 
obtained by generating 1,000 Monte-Carlo noise simulations 
for each line, in which Gaussian noise corresponding to the 
measured signal-to-noise ratio was added to the data and the 
fits repeated. The lower signal-to-noise present in the cores 
of the lines due to the fewer photons detected was taken into 
account during the Monte-Carlo process, under the assump- 
tion of Poisson statistics. Statistical uncertainties are given 
in table 3. 

In addition to these statistical uncertainties, there are a 
number of systematic uncertainties, primarily affecting the 
measured velocities. The small wavelength coverage of the 
UHRF (4-9 A) results in only a low number of arc lines for 
wavelength calibration which limits the absolute precision, 
particularly in the case of C2 for which only four arc lines 
were available. The magnitude of this effect can be esti- 
mated from the RMS of the wavelength fit (which used a 
cubic polynomial, quadratic for C2), which was always less 
than O.Olkms-i. Uncertainties in the atomic data listed in 
table 2 may affect the velocities e.g. the wavelength of the 
Call K line is known to be uncertain by about 0.02 km s-^ 
(Morton 2003). Velocity shifts between observations of a few 



tenths of a kms-^ have been seen in previous UHRF data 
(e.g. Crawford et al. 2000) and are thought to result from 
slight differences between the light paths of the telescope 
and arc lamp beams. The slight velocity offsets between the 
velocity components found in this work and those found by 
Crawford et al. (2000) are of similar magnitude. The quoted 
velocities are therefore systematically uncertain at around 
this level. Nevertheless, the multiple components for each 
species should be affected in an identical fashion by this sys- 
tematic uncertainty, so the relative velocities between com- 
ponents of a single species are unaffected, and so we do not 
incorporate this systematic uncertainty into table 3. 

Another potential source of systematic uncertainty is 
the definition of the continuum level. With the exception of 
Can (see section 4.2.2) there are no stellar counterparts to 
the interstellar lines. Around a thousand pixels of clear con- 
tinuum were present on both sides of every observed line. 
This greatly eases determination of the continuum level, 
which was fitted with a low-order polynomial (quadratic or 
cubic). The uncertainty in the continuum level is therefore 
neglible compared to the statistical uncertainties. 

Finally, we consider the systematic uncertainty intro- 
duced by the determination of the instrumental point-spread 
function (PSF). Inspection of the He-Ne laser and Th-Ar 
arc lamp spectra shows no evidence for a non-Gaussian 
PSF. The measurement of the width of the PSF (and 
hence instrumental resolution) carries an uncertainty of 
~ 0.001 km s- 1 , which translates to a systematic uncertainty 
of ~ 0.0006 kms-i on the measured h values. This uncer- 
tainty is negligibly small compared to the statistical uncer- 
tainties on 6, and affects all model components equally, so it 
has not been included in the uncertainties quoted in table 3. 

4.2 Atomic species 

Previous observers have modelled the main absorption 
around 9 kms- with either one or two separate com- 
ponents. Dunkin & Crawford (1999) and Crawford et al. 
(2000) found that their K I and Nal data required two com- 
ponents, whilst Crawford (2002) modelled only a single com- 
ponent in Cat and CH (but noted that these line profiles 
would also be consistent with the presence of two compo- 
nents, as found for the stronger lines). Whilst modelling the 
atomic species, it was found that a two-component model 
fitted to the narrower ~ 9kms-i feature was significantly 
better than a one-component model for all atomic species. 
This was confirmed through examination of the fit residu- 
als and via a statistical F-test, which found that the two 
component model was statistically better (in the sense of 
improvement in x^) at least the 95% confidence level. The 
similarities in the velocities and b values obtained for each 
species, which were derived independently (except for Kl, 
see below), reinforces the interpretation that a two compo- 
nent model is justified by the new high signal-to-noise data. 

For each atomic species, the main absorption around 
9kms-i contains a narrow component at ?; ~ 8.5kms-i 
with b ~ 0.7kms-i (hereafter component A) and a slightly 
higher velocity even narrower component at ~ 8.7kms-i 
with b ~ 0.3kms-i (component B). These components are 
consistent with the previously reported two-component b 
values of Crawford et al. (2000), though at slightly shifted 
velocity. All of the atomic species except Ca I also include a 
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Table 3. Results of the line profile analysis of the UHRF observations of interstellar absorption towards k Vel. The heliocentric velocity 
vq, Dopplcr broadening parameter b and logarithm of the column density A'^ (in cm~^) arc given for each model component. The total 
convolved line profiles are plotted in figures 1-5. T^' is the upper limit on the kinetic temperature calculated from equation 1 assuming 
no turbulent contribution to the b-value (see section 5.2). Heliocentric velocities are quoted with respect to the weighted means of the 
transitions for each species. Statistical uncertainties are quoted as Icr. 



b logAf t; 



(kms-^) (kms-^) (K) 

Cai Ag gQ4_|_Q Q33 0.851 ±0.049 9.412 ± 0.032 ITeO^l^T 

B8.762± 0.011 0.268 ±0.025 9.198 ±0.051 174l|^ 

Can -0.27 ±0.18 2.81 ± 0.21 10.123 ±0.033 

A8.500± 0.005 0.755 ±0.006 11.240 ±0.005 1381^22 

B8.802± 0.007 0.283 ±0.012 10.521 ± 0.024 194^^^ 

12.003 ±0.068 5.516 ±0.086 11.078 ±0.006 



Ki 



A8.200± 0.017 
B8.414± 0.008 

12.225^ 

13.70^ 



0.609 ± 0.036 
0.346 ±0.017 

1.31-' 

4.86-f 



10.141 ±0.014 
10.234 ± 0.009 

9.30 ±0.12 
9.979 ±0.050 



880: 
283 



-100 
+28 
27 



Nal 



^8.460 ± 0.007 
B8.720± 0.008 
12.225 ±0.039 
13.70 ±0.71 



0.694 ± 0.009 
0.286 ±0.018 
1.31 ±0.11 
4.86 ±0.55 



11.887 ±0.009 
11.327 ±0.032 
10.933 ± 0.053 
11.007 ±0.092 



671; 
114" 



-18 
-17 
-15 
-14 



CH 8.77 ±0.23 2.06 ± 0.56 11.470 ±0.051 

^ component A, ^ component B, 

f values held fixed during fitting (see text) 



broader absorption around ~ 12kms~^; this absorption was 
seen in Call and Kl by Crawford (2002), but only included 
in his Call model. The methodology adopted in modelling 
each species is described in more detail below. 



12.1 Cal 

There is no evidence for broad components in the Cal ab- 
sorption, so the model consists of the narrow components 
A & B only (see figure 1). Although the absorption only 
reaches a depth of 8% of the continuum, the high signal-to- 
noise data confirm the need for two narrow components. 



4.2.2 Can 

Calcium K is the only observed line which has a photo- 
spheric counterpart. The rotational velocity of k Vel is 
Vrot sini — 55kms^^ (Glebocki & Stawikowski 2000), so the 
stellar line is much broader than the interstellar features. 
Fortuitously, the orbital motion of the spectroscopic binary 
during our observations resulted in displacement of the stel- 
lar line to a velocity well clear of the interstellar absorption. 
In contrast to previous observations, the stellar line does not 
overlap the interstellar features in our data, which makes the 
continuum definition straightforward. The narrow absorp- 
tion around 9kms~^ clearly requires two components in the 
model (see figure 2). The broad higher- velocity absorption 
has been modelled with a single h = 5.516 ± 0.086 km s^^ 
component (Crawford 2002 also modelled this as a single 
component). There is no evidence for two broad components 
as found for Nal (see section 4.2.3 below). A new Call ab- 
sorption component at about Okms~^ is also found, which 



was regarded as part of the stellar continuum by Crawford 
(2002). 

4.2.3 Nai 

For this species it was found that two broad higher-velocity 
cloud components were required in addition to the two nar- 
row components A and B (see figure 3). These components 
were also found by Crawford et al. (2000), though the pa- 
rameters of the broad components differ in detail, typically 
at the 2-cr level (but less so when compared to the 1994 ob- 
servations). Our 12kms^^ component is slightly narrower 
and weaker, and our 14kms^^ component stronger and 
broader, than Crawford et al.'s corresponding components; 
however the total column density of the two broad compo- 
nents has in fact changed very little between 1994, 2000 and 
2006. Given their marginal significance, we therefore con- 
sider that the small differences between the modelled broad 
components might simply be due to the covariance between 
the fitted parameters, given the difficulty in defining accu- 
rately the velocities of such broad overlapping components. 
If the ~ 14kms~^ component is instead fixed to the veloc- 
ity found by Crawford et al. (2000), all other fitted parame- 
ters become quite consistent with the other epochs, for both 
broad components. There is therefore no evidence that these 
broad components are varying between epochs. 

There remains a slight discrepancy between the Nal 
model and data at the core of the narrow components, in- 
dicating that further substructure may be present. How- 
ever, models with additional components did not improve 
the fit, because the (superposed) blueward hyperfine com- 
ponent within the line cannot be reproduced without the 
redward hyperfine component becoming too deep relative 
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to the data. An investigation was undertaken to determine 
whether this behaviour could be explained by an error in 
the zero point of the spectrum, but adjusting the level of 
the zero point did not improve the fit. 

4.24 Ki 

For potassium a broad high-velocity absorption is visible 
(figure 4). This can be adequately modelled with a single 
component, but with a very high b value (~ 7kms~^). How- 
ever, because this absorption is seen to consist of two sep- 
arate components in Nal, which is expected to be similarly 
distributed to Kl (Pan et al. 2005), we have chosen to re- 
produce these two components in the K I model. Although 
the two components do not result in a statistically signif- 
icant improvement in the fit, the two-component model is 
more physically plausible, given the Nal model, and avoids 
depressing the continuum around the narrow components A 
and B. The broad components in the Kl model were fixed 
at the velocity and 6- values found for Nal, but were allowed 
to vary in column density. 

In an unconstrained fit, the narrow components (A and 
B) collapse to the same velocity, and the model becomes 
inconsistent with the velocity separation seen in the other 
atomic species. We therefore prefer a model for K I in which 
the components reproduce the velocity separation seen in 
earlier observations and in the other species. We iterated 
towards an acceptable model as follows; first the correla- 
tion between the distribution of Kl and Nal which was as- 
sumed for the broad components was extended to the nar- 
row components. The velocity separation between A and B 
was constrained to be consistent with both the other species 
and the two-component Kl model of Crawford et al. (2000). 
When the column densities and &-values in the ~ 9kms~^ 
component were fitted in this model, component A ap- 
pears to be weaker (log A'^a = 9.99) than the 1994 value re- 
ported by Crawford et al. (2000), with B somewhat stronger 
(logiVe = 10.32). Crawford et al. (2000) found that their re- 
sults were consistent with only the redward component (B) 
having varied between 1994 and 2000, so we also fixed the 
column density of component A to the 1994 value, to ob- 
tain the final model presented in table 3 (the uncertainties 
on the parameters of component A were obtained by fix- 
ing all other components at their fitted values, so are likely 
underestimated). In all unconstrained fits, component B is 
even stronger than tabulated here. However, the total col- 
umn density in both narrow components (A-l-B) is identical 
whichever approach is adopted, so that figure is robust even 
if the relative strengths remain uncertain. 

4-2.5 Comparison between atomic species 

Although the two main absorption components are broadly 
consistent between atomic species, their exact velocities, b- 
values and the velocity separation of the two components 
are not. These differences are most likely to be due to the 
systematic uncertainties discussed in section 4.1. Neverthe- 
less, the relative velocities between components of the same 
species are secure. 

The high- velocity components seen in Na I, K I and Ca 11 
have high b- values of up to 5 km s~^ . There are three possible 



explanations for such high 6- values: 

a) High temperature gas. We reject this possibility because 
the temperatures required (> 30,000 K) are inconsistent with 
the presence of substantial amounts of neutral Kl and Nal. 

b) Very turbulent gas, with Vt > 4kms~^. This would be an 
unusually high turbulent velocity, but cannot be excluded 
by the current data. 

c) Multiple blended absorption components, which cannot 
be distinguished at the current signal-to-noise ratio. This 
seems the most likely explanation (cf. Dunkin & Crawford 
1999). 

The differences in separation between the narrow com- 
ponents may indicate intrinsic unresolved structure within 
the absorbing material. Whilst the absorption has been 
modelled as two components, this does not necessarily im- 
ply that there are two physically separate clouds along the 
line of sight. In particular, the electron and number densities 
derived in section 5 are anomalously high for both compo- 
nents. There may exist a range of physical conditions within 
a single 'cloud', which are approximated by the two com- 
ponent model. The velocity separation between components 
A and B are larger for Can and Nal than for Cal (and 
Kl). This could potentially be explained if Call and Nal 
are present throughout the cloud, whilst Cal is present only 
in the denser regions and Kl is intermediate between Nal 
and Ca I (the distribution of these species in diffuse clouds 
has been discussed by Pan et al. 2005). If so, the physical 
parameters derived in section 5 are indicative of the range of 
conditions throughout the cloud, but will exclude the most 
extreme regions. This may also explain the imperfect fit of 
the model to the core of the Nal line. 



4.3 CH 

Whilst it might have been expected that CH would have 
a similar distribution to the atomic species, perhaps with 
enhancement in the cooler (and presumably denser) com- 
ponent B, no evidence for this has been found. There is no 
statistical justification for more than one component in the 
current data. The CH line profile was modelled with a single, 
broad (& = 2.06 kms~^) absorption component, albeit with 
a high uncertainty (ab = 0.56 kms~^), as shown in figure 5. 
The absorption profile has a central depth of less than 0.5% 
of the continuum. Although the velocity of the CH absorp- 
tion is closest to that of component B, the uncertainty and 
variations between species are large, so CH cannot be con- 
clusively identified with any component seen in the atomic 
species. We consider the breadth of the CH line profile in 
more detail in section 5.3. 



4.4 C2 

The continuum signal-to-noise ratio achieved in the region 
of the C2 lines was ~ 800 but no lines were found. To im- 
prove the chances of a detection, the expected positions of 
the five rotational lines were stacked in velocity space, but 
this also failed to reveal any absorption. A 3cr upper limit 
of 0.11mA was determined on each of the five lines. The 
corresponding upper limit on the total C2 column density is 
model dependent; we return to this point in section 5.8. 
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Table 4. Measured equivalent widths {W) of DIBs and atomic 
lines towards k Vel in January 1995 and June 2004. In the final 
column, 'Y' indicates a definite or probable stellar blend contam- 
inating the DIB, '-' indicates improbable or absent stellar blend 
and '?' is for indeterminate cases. Uncertainties and upper limits 
are CTcAA. 

W (mA) 



Jan 1995 June 2004 Blend? 



A5780 


10 ±4 


10 ±1 


7 


A5797 


< 1.5 


< 0.4 


7 


A5850 


< 1.2 


< 0.4 




A6196 


< 0.6 


1.1 ±0.2 




A6203 


< 6 


3±2 




A6283 




19 ±6 


Y 


A6613 


< 1.8 


3 ±0.5 




K I 7698 


t 


5.5 ±0.2 




Na I Di 


60.0 ±0.3 


62.2 ±0.2 





* Severe telluric contamination 

t Outside observed wavelength range 
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Figure 6. Spectra of the region around the A5780 and A5797 dif- 
fuse interstellar bands recorded towards k Vol. The upper trace is 
the observed spectrum in June 2004, the middle trace is the same 
region as observed in January 1995, and the lower trace shows 
residual intensities of the 2004 minus the 1995 data. The A5797 
DIB is not detected at either epoch, and there is no evidence for 
temporal variation of the A5780 DIB. 



4.5 Diffuse interstellar bands 

A search was performed for those DIBs in the observed wave- 
length regions with potentially detectable equivalent widths, 
by reference to the survey of diffuse interstellar bands by 
Hobbs et al. (2008). Five DIBs (AA5780, 6196, 6203, 6283 
and 6613) were detected towards k, Vel in the 2004 data; 
only A5780 was detected in the lower signal-to-noise 1995 
data (see figures 6-8). As far as the authors are aware, these 
are the first reported detections of DIBs towards this star, 
so there are no other epochs for comparison. 

The measured DIB equivalent widths are given in ta- 
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^ 1.005 
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Figure 7. As in figure 6, but for the A6196 and A6203 DIBs. Both 
DIBs are detected in 2004 only. 
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Figure 8. As in figure 6, but for the A6613 DIB. The line at 
6610 A is due to stellar Nil, whilst the DIB is detected in 2004 
only. 



ble 4, along with Nal Di which was present in both obser- 
vations, and the Kl 7698 line which was within the wave- 
length coverage in 2004 only. Uncertainties and upper limits 
were calculated as UcAA, where (Jc is the continuum noise 
level and AA the DIB FWHM. This simple method is fairly 
conservative and probably over-estimates the uncertainties, 
particularly for the broadest bands. The A6283 DIB over- 
laps the telluric O2 a-band, which could only be adequately 
removed from the 2004 data. Blends with stellar lines were 
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Figure 9. Comparison of observed equivalent widths of interstel- 
lar absorbers towards re Vel at different epochs. Equivalent widths 
of the interstellar absorption in 1994, 2000 (Crawford et al. 2000), 
2002 (Crawford 2002), 1995, 2004 and 2006 (this work) are plot- 
ted against stellar displacement on the bottom axis, and epoch in 
years on the top axis. Data have been normalised with respect to 
the 2006 values and vertical offsets of 0.2 between species applied 
for display; lines have been added to guide the eye. 
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Figure 10. As figure 9, but for the combined column density of 
the interstellar absorption at ~9kms~^. Only the UHRF data 
from 1994, 2000, 2002 and 2006 are plotted; error bars for 1994 
and 2000 are probably overestimates by factors of 5-10 (sec sec- 
tion 5.1.2). 



determined by reference to the stars v Sco and a Pav (both 
B2 IV) , which were observed on the same night as the 2004 
data. For A6283, a narrow stellar line at 6287 A was removed 
via a Gaussian fit; there are possible weak stellar blends with 
A5780 and A5797. The displacement in the continuum adja- 
cent to 5780 is likely to have been introduced as a result of 
weak stellar features in the telluric standards. 

As a comparison between the two epochs, the spectra 
around A5780 (the only DIB detected at both epochs, and 
including the region where A5797 would be expected) were 
subtracted and the residuals evaluated. This procedure is 
shown in figure 6; within the noise there has apparently 
been no change in A5780 between January 1995 and June 
2004. From the residuals and the uncertainties in the mea- 
sured equivalent widths, we determine an upper limit on the 
change in A5780 between the two epochs of <40%. 



5 DISCUSSION 

5.1 Comparison of observations at different 
epochs 

5.1.1 Equivalent widths 

A model-independent comparison of the measured equiva- 
lent widths of each atomic and diatomic line from all avail- 
able epochs since 1994 is presented in figure 9. For compari- 
son with the earlier measurements from Dunkin & Crawford 
(1999), Crawford et al. (2000) and Crawford (2002), we have 
consistently compared equivalent width measurements made 
over the same velocity range. 

For Kl, the combined equivalent width of only the nar- 
rower ~ 9kms^^ components has been plotted, to be con- 
sistent with the earlier data; this was measured as 5.09 ± 
0.07 mA in the 2006 data by first dividing out the model fit 
to the broad, higher velocity Kl components. These broad 
components are unavoidably blended with the narrow com- 



ponents in the lower-resolution 2004 UCLES data, leading 
to that datum appearing too high in figure 9. To check this, 
the 2006 Kl observation (including the broad components) 
was degraded to the same resolution as the 2004 observa- 
tion, and the equivalent width remeasured. The resulting 
value was 5.9 mA, which is somewhat lower than the total 
7.49 mA measured from the UHRF data (table 1) but still 
higher than the 5.09 mA of the narrow components alone. 
The equivalent width measured at the degraded resolution 
is thus similar to that measured in the 2004 UCLES data 
(see table 4); so the apparently inflated 2004 datum could 
be an effect of the lower resolution. 

There are significant increases in equivalent width for 
the K I and Ca I lines. W{K i) increased by 82 ± 7% between 
1994 and 2006, whilst VF(Cal) increased by 26±9% between 
2002 and 2006. The equivalent width of Nal barely changed 
over the same time periods (4 ±4% between 1994 and 2006, 
and 2 ± 4% between 2002 and 2006), ahhough this may 
be due to the near-saturation of the Nal absorption, which 
lies on the fiat part of the curve of growth. The measured 
Nal equivalent width also includes the broad, higher velocity 
components which are not suspected to be variable, given 
the lack of any significant differences between the models 
of the 1994, 2000 and 2006 data (see section 4.2.3). Call 
shows no significant increase (6 ± 4%) in equivalent width 
between 1994 and 2006, at least when matching velocity 
components are compared; the 2006 measurement plotted 
here (20.6 ± 0.2 mA) excludes the ~ Okms~^ component, 
but the (apparently low) 2002 datum is almost certainly 
compromised by interference from stellar Call K, which was 
well separated from the interstellar Ca ll K line in the 1994 
and 2006 data. The increase in CH equivalent width is not 
significant at 19 ± 22% between 2002 and 2006, but nor are 
the data inconsistent with CH matching the increases seen 
in K I and Ca I. 
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5.1.2 Column densities 

The problems caused by the near-saturation of the Nal 
line and the inclusion of the broad ~12kms^^ component 
may be avoided by instead comparing the modelled column 
density at each epoch, and restricting comparison to the 
~9kms~^ component (s) only. This is impractical for the 
Mt. Stromlo and UCLES Kl data, in which the narrow and 
broad components are not resolved, so they are not consid- 
ered further. The individual (A and B) components in the 
UHRF data cannot always be directly compared with the 
previously published column densities because the absorp- 
tion has been modelled with different numbers of compo- 
nents (e.g. for Cal and Can). However, all the remaining 
observations were made with the same instrument at essen- 
tially the same spectral resolution (850,000-900,000), and 
the total column densities (A + B) in the 9kms~^ compo- 
nents are robustly determined. 

In figure 10, the components reported in Crawford et al. 
(2000); Crawford (2002) and in this work are compared by 
summing the column densities of the one or two components 
which make up the absorption at ~9kms^^. Uncertainties 
in the total column density obtained simply by summing 
individual components in quadrature are overestimates, be- 
cause in overlapping lines the column densities of the two 
components are not independent quantities. The 2002 obser- 
vations are not affected, because they only utilised one com- 
ponent. For the 2006 observations, the uncertainties were ob- 
tained by interrogating the Monte-Carlo error simulations, 
determining the sum of column densities at each iteration, 
and taking the 68% range. This process revealed a strong 
anticorrelation between the components, and the total col- 
umn densities were found to be logA'AB(Cal) = 9.62±0.01, 
logiVAB(Can) = 11.316±0.002, logiVAB(Ki) = 10.47±0.01 
and logiVAB(Nai) = 11.985 ± 0.003. For the 1994 and 2000 
observations, only summing in quadrature is possible; by 
comparison with the 2006 observations these error bars are 
likely to be overestimates by factors of 5-10. 

The differing treatments of the higher velocity absorp- 
tion may still affect this comparison. For example, the broad 
Kl components partially overlap the ~ 9kms~^ absorption, 
but were not included in the model of Crawford (2002); this 
has the effect of slightly enhancing the apparent K I column 
density in 2002. 

Figure 10 shows that the increase in Kl column den- 
sity first noted by Crawford et al. (2000) has continued un- 
til 2006: iV(Kl) has increased by 82tJ°% between 1994 and 
2006. For the first time, Cal is found to have increased - by 
32 ± 5% between 2002 and 2006 (a slightly higher rate than 
Kl over the same period). Can and CH show no statistically 
significant changes between 2002 and 2006, and nor does Nal 
over the period 1994-2006. Cal and Call are therefore seen 
to vary at different rates to each other, indicating a change 
in the ionisation conditions. 
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Figure 11. Observed and simulated CH line profiles. Upper his- 
togram is the observed CH profile as in figure 5, middle histogram 
is the 'thermal' model, and the lower histogram is the 'turbu- 
lent' model, assuming CH shares the same distribution as Kl 
(see section 5.3 for details of the models). The spectra have been 
vertically ofi'set by 0.01 for display; the dashed lines illustrate the 
continuum levels. The observed line profile is significantly broader 
than the turbulent model, and the A-doublet is not visible in the 
observed spectrum. 

mechanisms which act to broaden interstellar lines are ther- 
mal Doppler broadening and turbulent bulk motions of the 
gas. These are related to the Doppler broadening parameter 
b through: 

where fes is the Boltzmann constant, Tk is the kinetic tem- 
perature, m is the atomic (or molecular) mass of the species 
and vt is the RMS turbulent velocity. Taking the limiting 
case of zero turbulent velocity, rigorous upper limits on the 
kinetic temperature T"' may be obtained for each species, 
and are given in table 3 when T^"' < 2,000 K. 

In principle it should be possible to determine Tk and 
Vt separately by solving equation 1 for species with differing 
values of m. However, there are no systematic differences 
in b values between species in component B, over a range 
m = 23 (Nal) to m = 40 (Cal and Call). In component 
A, statistically significant differences in b are found between 
species, but with Kl the narrowest and Cal the broadest. 
This is difficult to reconcile with equation 1, which implies 
that the least massive species (Na) should have the highest 
6 value. The most likely explanation for these discrepancies 
is that the species are not entirely physically coincident and 
reside within different parts of the cloud, perhaps in regions 
with differing vt. If Vt makes a significant contribution to the 
6 values, then the kinetic temperature could be substantially 
lower than the derived upper limits. 



5.2 Line widths and temperature 

Although one of the aims of the UHRF observations was to 
measure the kinetic temperature from the rotational pop- 
ulations of C2, the lack of a detection of this molecule ne- 
cessitates the use of another method. In addition to the in- 
strumental resolution and the natural line width, the two 



5.3 CH distribution 

Visual inspection of the CH spectrum shows no sign of the 
expected A-doublet line profile, which causes the relatively 
high value of b in the single-component vapid model. In 
order to determine if CH has the same line profile as the 
atomic species but the A-doublet is being hidden by the 
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noise, simulated CH line profiles were constructed. Two com- 
ponents were assumed, with velocity separation and col- 
umn density ratio equal to that observed in Kl, which is 
known to be distributed similarly to CH in diffuse clouds 
(e.g. Welty & Hobbs 2001; Pan et al. 2005), and a total col- 
umn density equal to the CH model. Two limiting cases 
for b were considered: vt dominant (the 'turbulent' model, 
Tk = 0), and dominant (the 'thermal' model, vt =0). For 
the turbulent model we assume b values equal to those of K I 
(table 3). Because CH (m = 13) is substantially less massive 
than K (m = 39), the thermal model adopts b values of 1.06 
and 0.60 km s~^ i.e. a factor of y/3 (= y/39/13) higher than 
the turbulent model. These simulated line profiles were then 
degraded to the obtained noise level using clean sections of 
the observed continuum near the CH line, which avoids the 
need to make assumptions about the noise statistics. 

The simulated spectra are compared with the observa- 
tion in figure 11. The A-doublet is clearly visible in the tur- 
bulent model, but barely discernible in the thermal model. 
A simulation in which only component A is 'turbulent' is not 
very different from the combined thermal model illustrated: 
it is the degree of turbulent broadening in the narrower com- 
ponent B which is primarily constrained by this analysis. It 
is therefore plausible that the CH shares a similar distribu- 
tion to Kl but the line profile is being masked by noise. If 
so, the turbulent contribution to the CH line profile is low, 
at least in component B. This suggests that the temperature 
of component B cannot be much lower than the upper limit 
derived in section 5.2. 

We also constructed simulated spectra based on the Ca I 
distribution, which is more heavily concentrated in compo- 
nent B than K I. In those simulations the A-doublet is clearly 
visible even in the thermal model, which is inconsistent with 
the observations. This supports the assumption made above 
that CH is distributed more similarly to Kl than Cal. 



5.4 Cal / Call ratio and electron density 

At the low densities prevalent in the diffuse ISM, the balance 
between the various ionisation states of a single atom is set 
by the competition between photoionisation which promotes 
the atom to a higher ionisation state, and recombination 
with an electron which demotes it to a lower one. For the 
neutral and first ionised state in equilibrium, this may be 
expressed as: 

r(Xi)n(Xi) = a^(Xi,T)n(Xn)ne (2) 

where P is the photoionisation rate, ar the radiative recom- 
bination rate coefficient, the electron number density and 
n(X) the number density of species X. As observations pro- 
vide the column density A'' and not the number density n, it 
is necessary to assume a constant density along the line of 
sight, allowing n to be replaced by A*'. Under this assump- 
tion, and for the specific case of calcium, equation 2 may be 
rewritten as: 

_ P(Cal) jV(Cal) 
" Q,(Cal,r) Ar(Call)' ^ ' 

Calcium is an unusual case because both its first and 
second ionisation energies (6.11, 11.87 eV respectively) are 
lower than the Rydberg energy (13.6 eV) , so in diffuse clouds 



ionisation to Cam is possible. There is then a second ionisa- 
tion equilibrium between Can and Cam, though this does 
not invalidate the electron densities derived from equation 3. 
Welty et al. (2003) calculated the relative populations of 
Cal, Can and Cam and found that the dominant ionisa- 
tion state is Can in any situation where Cal is abundant 
enough to be detected. 

Equation 3 is an approximation which assumes that no 
other processes cause significant conversion between ionisa- 
tion stages. Plausible processes include cosmic-ray ionisa- 
tion and charge exchange with dust grains or PAHs (e.g. 
Liszt 2003). Welty et al. (2003) have discussed this issue 
in the context of calcium ionisation ratios, and concluded 
that equation 3 remains valid if the fractional ionisation 
rie/nn < 2 x 10~*, which is the level provided by com- 
plete ionisation of gas-phase C I to C n (assuming standard 
depletions, see also section 5.7). Those authors also conclude 
that the assumption of photoionisation is most valid when 
He is high. Since the observations towards k Vel indicate a 
cloud which is much cooler and/or more dense than usual 
for the diffuse ISM, and as the low extinction present in the 
cloud'' favours photoionisation, the assumptions leading to 
equation 3 are expected to be valid in this case. 

The electron density may then be determined from the 
column densities, given the value P/or = 66 for Ca I given by 
Welty et al. (2003, their table 3)^. The electron densities are 
then found to be Tie, a ~ 0.98 ± 0.12 cm~'^and ne,B ~ 3.14 ± 
0.77 cm^^ for components A and B respectively. These are 
very high values for diffuse clouds. For comparison, Crawford 
(2002) found rie — 0.97 cm~^ towards k Vel using a single- 
component model. 

Some observers have found that electron densities de- 
rived from calcium are sometimes higher than those derived 
from other species along the same lines of sight. Welty et al. 
(1999) used ionisation ratios of nine different elements to 
determine Ue towards 23 Ori, with calcium resulting in the 
highest estimate, whilst Sonnentrucker et al. (2003) found 
that calcium gave the highest rie (out of four measured ioni- 
sation ratios) towards HD 185418. This effect has only been 
reported towards a handful of stars, and the cause is un- 
known (see Welty et al. 2003 for a discussion) . There is some 
circumstantial evidence to support the electron density de- 
rived from calcium: the changes seen in Kl are about the 
same rate as those seen in Cal (section 5.1), suggesting that 
both are responding to the same change in (observed) elec- 
tron density; and the total number densities derived inde- 
pendently from the electron density and the chemical model 
(sections 5.7 and 5.9) are in agreement with each other. 

5.5 Cal / Kl ratio and depletions 

Whilst the Cal/Can ratio provides a measure of the elec- 
tron density (under the assumption of constant tempera- 
ture), some line ratios are sensitive to other properties of the 

^ -B(a_v) = 0.10, so Ay = 0.16 halfway through the interstellar 
material along the line-of sight, assuming Ry =3.1. The trans- 
verse extinction is unknown but probably even lower. 
^ This value assumes T = 100 K. For calcium, the dependence of 
ctr on T is o oc (T/lOO K)-0 S83 (Pequignot & Aldrovandi 1986). 
r is independent of T. 
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interstellar gas. Because the depletion of potassium into dust 
grains is thought to be fairly constant, the Ca l/K I ratio may 
be used to estimate the depletion of calcium, which is typ- 
ically very high at the low temperatures in the diffuse ISM 
due to its high condensation temperature [Tc = 1,659K, 
Lodders 2003). Welty et al. (2003) found that Ar(Cai) and 
A'^(K l) are correlated, with the slope of that correlation aris- 
ing from the difference in depletion behaviours. Our mea- 
sured ratios fall very close to the mean relation determined 
by Welty et al. (2003), suggesting that there is no unusual 
enhancement of Cal along this line of sight. Adopting the 
same assumptions as for the calcium ionisation balance, and 
assuming that Call is dominant, the relative populations of 
Cal and Kl are given by: 



iV(Cai) ^ (r/gQcai 5{Ca.)ac. 
iV(Kl) (r/a^Ki 5(K)aK 



(4) 



where 5{X) is the gas-phase fraction and ax the cosmic 
abundance of species X. Taking the standard abundances, 
5(K) = 0.2, and photoionisation and recombination rates 
given by Welty et al. (2003) this becomes 



^•(Cal) 
iV(Ki) 



12.4(5(Ca). 



(5) 



This gives estimates of the calcium gas-phase fraction in 
components A and B of 5(Ca)A = (1.51 ± 0.23) x 10"^ and 
5(Ca)B = (7.42 ± 1.45) x 10-^ which may be equivalently 
expressed as D(Ca)A = -1.82±0.13 and L>(Ca)B = -2.13± 
0.17, where D are the logarithmic depletions. Cardelli et al. 
(1991) investigated the depletion of calcium in diffuse clouds, 
and concluded that it is very sensitive to the local density, 
increasing in denser regions. This is in accordance with our 
result of a higher calcium depletion in component B (which 
is found to be approximately a factor of two denser than 
component A in section 5.7). 



5.6 Molecular hydrogen fraction 

The Cal/Kl ratio is also correlated with the fraction of 
hydrogen in the form of H2, though the relation is not 
particularly tight (Welty et al. 2003). Previous observa- 
tions have only resulted in upper limits for the Hi and 
H2 column densities towards k Vel (Jenkins 2009, whose 
results were based on the observations of Bohlin et al. 
1983 corrected for stellar contamination). The measured 
ratios are log(Ar(Cai)/Ar(K i))a = -0.729 ± 0.035 and 
log(Ar(Cal)/Ar(Kl))B = -1.036 ± 0.052. These correspond 
to upper limits for the fraction of hydrogen in the form 
of H2 of /(H2) < 10"^ in both components (Welty et al. 
2003, their figure 4). These are very low molecular hydro- 
gen fractions for the diffuse ISM, especially for regions with 
enhanced density. 

A more direct estimate of the molecular hydrogen frac- 
tion could in principle be obtained from the standard rela- 
tions between E(^b-v) s^nd Hi, and between A''(CH) and H2. 
Adopting E(^B-v) ~ 0.10 (Cha et al. 2000) with an uncer- 
tainty^ of 0.02, equation 7 of Burstein & Heiles (1978) pre- 
dicts log Ar(Hi) = 20.9 ± 0.2. Taking log(A/'(CH)/iV(H2)) = 



-7.46 ± 0.21 from Sheffer et al. (2008)**, we obtain an esti- 
mate of logiV(H2) = 18.93 ± 0.22. This implies a molecular 
fraction of /(H2) = ,^if^'(H, ) = (2.1±1.6) x 10-^ which 
is in line with the previous estimate. 

Both of these estimated column densities are more than 
one dex higher than the upper limits from the UV obser- 
vations (Jenkins 2009). We also note that Welty & Hobbs 
(2001) presented empirical relations between Kl, Nal and 
total hydrogen column density for typical diffuse ISM lines 
of sight, but for k Vel these would again imply at least an 
order of magnitude more hydrogen than the observed upper 
limits. Apparently the standard relations are not valid for 
the line-of-sight towards k Vel. It is possible that the column 
densities have changed by this amount between the epochs 
of the UV observations (late 1970s, Bohlin et al. 1983) and 
the optical observations (the E(^b-v) measurement is based 
on Tycho photometry so has an effective epoch of 1991.5, 
whilst the spectroscopic measurements are from 2006). Al- 
ternatively the ISM towards «; Vel could be out of chemical 
equilibrium (Bell et al. 2005) and thus these standard rela- 
tions might not apply. 



5.7 Physical dimensions and density 

The proper motion of k Vel is 15.5 mas yr^^ (Ferryman 
1997), so at a distance of 165 pc the transverse velocity is 
2.56 AU yr~^ (12.1 km s~^). This gives a stellar displacement 
of 10.1 AU between the 2002 and 2006 observations, and 
29.7 AU between 1994 and 2006 (shown on the axes of fig- 
ures 9 and 10). The size of the binary orbit has been deter- 
mined to be ai sini = 0.48 AU, where a\ is the semi-major 
axis of the orbit of the primary, and i the orbital inclination 
(Pourbaix et al. 2004). Since neither the inclination nor po- 
sition angle of the orbital ellipse have been determined, it 
is impossible to include the effect of this on the transverse 
displacement, but the error thus introduced is likely to be 
rather small. 

The displacement of the line-of-sight at the distance 
of the star provides an upper limit on the displace- 
ment at the distance of the foreground absorbing cloud. 
Dunkin & Crawford (1999) have argued that the cloud is 
likely to be located close to the star, at least 150 pc from 
the Sun, based on comparison with CO maps and optical ab- 
sorption along other nearby lines-of-sight (see also Crawford 
1991). This would place the cloud just beyond the edge of 
the Local Bubble (Welsh et al. 2010). Any transverse mo- 
tion of the interstellar cloud would modify the displacement 
of the line-of-sight through the cloud by the vector sum of 
the motions of the star and cloud. It is impossible to de- 
termine the transverse velocity of the cloud, but reasonable 
assumptions suggest that the derived cloud size is accurate 
to within a factor of 2-3, in either direction. 

If the electron density is dominated by electrons re- 
leased through the photoionisation of atomic carbon (the 
assumption adopted by Crawford 2002), the electron den- 
sity allows an estimate of the total number density to be 
obtained. Assuming a carbon depletion of 60% implies a gas- 
phase carbon abundance of 1.4 x 10~*nH (Sofia et al. 1997), 



^ The uncertainty on -E(g_v) is taken to be equivalent to a mis- 
elassification by one spectral subtype. 



* Note that this requires extrapolating the Sheffer et al. (2008) 
relation to lower column densities than they actually observed. 
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which leads to estimates of tia > 7 x 10"^ cm"'' and ne > 
2x 10'' cm~"^. These are orders of magnitude higher than typ- 
ical values for diffuse or translucent clouds (Snow & McCall 
2006). They are also higher than the densities inferred by 
some small-scale structure studies (e.g. Pan et al. 2001), but 
those studies relied upon simple chemical models to deter- 
mine the density. Bell et al. (2005) demonstrated that the 
chemistry can be time-dependent in small-scale structures, 
whilst the atomic species react much more rapidly to chang- 
ing conditions. 

Welty (2007) combined optical and UV observations 
(including C I and C ll) of the small-scale structure towards 
HD 219188. It was inferred that the changing column den- 
sities reflected enhanced recombination due to a raised elec- 
tron density, the source of which was the ionisation of hy- 
drogen. If this effect is also present in the line-of-sight to- 
wards K Vel, the densities derived above will be overesti- 
mates. Since hydrogen is thousands of times more abundant 
than carbon, even a small fractional ionisation of hydrogen 
could reduce the density estimate by an order of magnitude 
or more. An ionisation fraction of > 1% would be required 
to reduce the density to levels typical of the diffuse ISM. A 
physical cause of the ionisation would also be required, such 
as intense cosmic ray bombardment or penetration of soft X- 
rays through the relatively low-extinction gas (Tielens 2005), 
though the latter might plausibly arise from the hot gas in 
the Local Bubble (Frisch et al. 2011). 

An estimate of the hydrogen column density may be 
derived from the observed metal column densities and de- 
pletion estimates. Taking log[Ca/H] -I- 12 = 6.34 from 
Asplund et al. (2009), D(Ca) estimates from section 5.5 
and the observed Cal -I- Call column densities, we obtain 
logArA(H) ^ 18.7 and logAfB(H) ^ 18.3, so logiVAB(H) ^ 
18.9. This is consistent with the observational upper limits 
of logiV(Hl) < 19.5 and logiV(H2) < 17.7 (Jenkins 2009). 
Taking these column density estimates and assuming a den- 
sity of 10* cm~^ gives an estimated total path length of 
~ 50 AU. This length estimate rests upon a series of as- 
sumptions, in particular the derived calcium depletions and 
neglecting Cam, so is probably only indicative to within an 
order of magnitude. 

5.8 C2 column density 

The relative populations of the low- J rotational states of C2 
depend upon the local density and kinetic temperature. In 
order to determine an upper limit on the total C2 column 
density it is necessary to combine the limits on the equiv- 
alent width for lines with differing J derived in section 4.4 
with the density and temperature determinations from sec- 
tions 5.2 and 5.7. Since C2 is expected to preferentially form 
in denser regions, we adopt n = 2 x 10* cm""' and assume 
Tfe = 100 K. 

The relative populations of each J level were calculated 
using B. McCall's online calculator'^ which uses the method 
of van Dishoeck & Black (1982). The populations of each 
level were combined with the oscillator strengths from ta- 
ble 2 and the calculated total population (i.e. integrating 
over all J) to determine the maximum total column density 

^ http://dib.uiuc.edu/c2 



which is compatible with the equivalent width upper limits. 
The resulting 2a upper limit is iV(C2) < 5.6 x lO" cm"^. 

Adopting a higher density or lower temperature would 
result in a somewhat tighter upper limit. At — 100 K 
the calculation is rather insensitive to n: if n changes by an 
order of magnitude the calculated upper limit differs by less 
than 5%. Varying Tk by 20 K changes the upper limit by 
~ 10%. 

5.9 Comparison with chemical models 

The chemical models of Bell et al. (2005) made predictions 
for the column densities of CH, C2, Cal and Call (plus 
OH and CO) along a line of sight through an SSS filament. 
The physical parameters were based on those derived by 
Crawford (2002) for k Vel, namely an extended structure of 
length 100 AU along the line-of-sight and transverse depth 
of 10 AU. Results were presented for a grid of cloud densities 
and times. 

Because CH shows no evidence for more than one cloud 
component and C2 has only an upper limit over the to- 
tal line-of-sight, for comparison with the chemical model 
column densities were initially summed over both compo- 
nents A and B. The total column densities adopted in this 
case"' were iVAB(CH) = 2.95 x lO" cm^^ iVAB(C2) < 
5.6 X 10" cm"^ (2(7), iVAB(Cal) = 4.17 x lOVm'^ and 
iVAB(Call) = 2.07x10" cm"^ (cf. section 5.1). With the ex- 
ception of Call, these column densities are fully consistent 
with the Bell et al. results for n = 5 x 10*cm~^, t = 10- 
50 yr. Models with n — 10^ cm~^ are excluded by the C2 
upper limit, whilst those with lower density are excluded by 
the Ca I column density. The observed Ca 11 column density 
is moderately higher than the model, by a factor of ~ 4, but 
this is within the uncertainties of the models. 

In an attempt to include material only in the densest 
regions, we have also compared the Bell et al. (2005) mod- 
els to component B alone - this component was suspected 
to be variable by Crawford et al. (2000). It then becomes 
necessary to assume that all of the CH arises within com- 
ponent B (but see the discussion in section 4.3), and the 
high column density of Cal in component A cannot be ex- 
plained. Bearing in mind these caveats, the CH and C2 col- 
umn densities remain as above whilst the calcium column 
densities are A^B(Cal) = 1.58 X 10^ cm"^ and iVB(Call) = 
3.32 X 10^" cm~^. All the observed column densities are then 
consistent with the Bell et al. model for n = 5 x 10* cm"'', 
t = lOyr. 

The chemical models thus provide confirmation of the 
high number density in the cloud towards k Vel, and are 
most consistent with densities only slightly higher than the 
lower limits derived in section 5.7. The predicted C2 column 
density in the most consistent model is roughly a factor of 
three times lower than our observed upper limit. The models 
determined the local temperature through a self-consistent 
treatment of heating and cooling processes. The range of 
temperatures in the models was ~ 30-115 K, which are con- 
sistent with the upper limits derived in section 5.2. It is dif- 

Although components A and B were not observed in CH and 
C2, we assume that these molecules are confined to those cloud 
components. 
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ficult to interpret the time dependent aspect of the models, 
because the structure cannot have formed instantaneously. 
Nevertheless, the models predict that for n = 5 x 10* cm^'^ 
CH and Ca ll would be essentially constant over time whilst 
Cai would increase, as observed (cf. figure 10). 

Bell et al. set n(H2)/n(Hl) = 0.4, which is more than 
an order of magnitude higher than the molecular hydro- 
gen fraction inferred in section 5.6. Low molecular hydrogen 
fractions would act to reduce the column densities of the 
molecular species because the formation pathways all in- 
clude reactions with H2. The estimates in section 5.6 might 
therefore be misleading, as the correlations may not hold 
in such high-density regions. Alternatively, the local den- 
sity and/or H2 fraction could be highly variable along the 
line-of-sight. A large number of such variations in local den- 
sity (and to an extent temperature) have been proposed by 
Cecchi-Pestellini et al. (2009) as a solution to several prob- 
lems in diffuse cloud chemistry. If this were the case, it would 
be possible for CH to be formed in regions of relatively high 
H2 density, but for the atomic lines to trace less-dense ma- 
terial elsewhere, therefore reducing the total molecular frac- 
tion along the line of sight. In the absence of direct ob- 
servations, this contradiction in ,f (H2) cannot currently be 
resolved. 



5.10 Diffuse interstellar bands 

The measured DIB equivalent widths are far lower than ex- 
pected given the reddening, and only a few DIBs have been 
detected. Based on 133 lines of sight, Friedman et al. (2011) 
found a mean 

W{\h7iO)/E(B-v) value of 505 mA mag 
However, towards n Vel the reddening is E(^b-v) = 0.10 ± 
0.02 but W(A5780) is 10 ± 1 mA, i.e. five times lower than 
would be typical for the Milky Way ISM. The relative weak- 
ness of A5780 indicates a lower abundance of its carrier in 
the dense material towards k Vel than in the ambient ISM, 
which may indicate a density dependence of the (unknown) 
carrier formation mechanism. 

Those DIBs which have been detected towards k Vel 
are members of Krelowski & Walker (1987) group 2, whilst 
the DIBs of group 3 are absent. The equivalent width ra- 
tio of A5780 (a group 2 DIB) to A5797 (group 3) is at 
least W(A5780)/M^(A5797) > 22 (la), which is amongst 
the highest known: typical values of this ratio are 1-4. 
As far as the authors are aware, the highest ratios previ- 
ously reported in the literature are W (A5780)/W(A5797) > 
9 (Ict) towards the star eCassiopeiae by Friedman et al. 
(2011), and Vy(A5780)/W(A5797) = 8.2 (no uncer- 
tainty given) towards HD 20336 by Galazutdinov et al. 
(1998)". Krelowski et al. (1999) found a correlation be- 
tween H/(A5797)/Vy(A5780) and W iCB) / E^g-v)- Although 
the value of W{CR)/E(^b~v) = 2.5 ± 0.6 mA mag ^ mea- 
sured towards k Vel is much lower than any measured by 
Krelowski et al., the upper limit on (A5797)/W(A5780) 
does appear to fall on or near this relation. 

The observed A5780/A5797 ratio corresponds to a cloud 

Porcoddu ct al. (1992) report five sightlines with high A5780/ 
A5797 ratios, and Galazutdinov et al. (1998) three more, but nei- 
ther set of authors report the value of the ratio nor provide nu- 
merical upper limits on VF(A5797). 



of extreme a-type (Krelowski et al. 1992), which is generally 
interpreted as a low density, high UV radiation field environ- 
ment. The additional DIBs and diatomic molecules present 
in <^-type clouds (generally interpreted as cooler, denser and 
less UV-irradiated than a-type) are extremely weak or ab- 
sent. However, the analysis above indicates that the inter- 
stellar cloud(s) toward(s) n Vel is/are relatively dense but 
unshielded from UV radiation. The weakness of A5780 with 
respect to E(^b-v) and the extreme cr-type line-of-sight sug- 
gests that unusual conditions are present in the cloud to- 
wards K Vel. Cami et al. (1997) found similar properties 
(weak DIBs per unit -E(_b-v), extreme A5780 to A5797 ratio, 
weak or absent diatomics) towards targets in Orion, where 
the conditions are also fairly dense but highly UV irradi- 
ated. Snow et al. (1995) found systematically weak DIBs in 
reflection nebulae, where conditions are even more extreme 
in terms of high UV field and high electron density. 

It has been suggested that the carrier of the A5797 
DIB is a neutral species which is easily ionised, whilst 
the A5780 carrier does not easily change ionisation state 
due to being either an already ionised species or a neutral 
which is highly resistant to ionisation (e.g. Cami et al. 1997; 
Krelowski et al. 1998) Because the A5780 intensity over a 
large sample correlates very well with A*'(H l) it has also been 
suggested that the carrier has an ionisation energy similar 
to that of hydrogen (Friedman et al. 2011). Whilst the hy- 
pothesis on the nature of the A5797 carrier is in accordance 
with the widely observed enhancement of A5797 in less UV- 
irradiated environments, it appears to be at odds with the 
lack of A5797 detection towards k Vel given the high electron 
densities derived in section 5.4. 

High electron densities should lead to enhanced recom- 
bination of molecular cations to form neutral molecules, 
as well as the observed neutral atomic species. But this 
would result in relatively strong A5797, opposite to what 
is observed towards re Vel. Electronic recombination with 
a molecular cation does not necessarily produce the neu- 
tral through radiative association because the molecule may 
lose energy through fragmentation instead of emitting a pho- 
ton. But if such fragmentation pathways are significant for 
the carrier of A5797 this should result in destruction of the 
carrier in environments with high electron density; this is 
incompatable with the hypothesised enhancement of A5797 
in dense environments due to enhanced recombination of a 
cation. 

It appears that the carrier of the A5797 DIB (and by 
extension, the other DIBs in Krelowski & Walker group 3) 
is more sensitive to UV irradiation than density or degree of 
ionisation. The observations of re Vel may indicate that the 
carrier of A5797 is dissociated by UV radiation more rapidly 
than the carrier of A5780. 

Although the Kl line was not observed in 1995, by in- 
terpolating figure 9 it is inferred that the equivalent width 
of the Kl line increased by « 80% between the 1995 and 
2004 observations; the limit on the change in A5780 DIB is 
< 40%. If the carrier of A5780 were to have increased in 
proportion with the K I line, a change would have been seen 
between 1995 and 2004, albeit with low statistical signifi- 



Snow et al. (1995) argued against already ionised species, in- 
stead favouring neutral molecules. 
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cance. Cordiner et al. (2006) and Cordiner et al. (in prep.) 
found that DIB carriers (including A5780) trace the same 
SSS as atomic lines such as Kl. The fairly loose upper limit 
on the variation of A5780 towards k Vel is insufficient to 
draw firm conclusions on this point. 



6 CONCLUSIONS 

New ultra-high resolution observations of interstellar ab- 
sorption towards k Vel by the species K I, Na I, Ca I, Ca ll and 
CH have been obtained. A search for C2 was undertaken, but 
no lines were detected. Best-fitting model line profiles were 
computed and two main narrow (b < 0.8kms~^) absorption 
components were identified. Comparisons with observations 
taken at earlier epochs show increases in interstellar K I and 
Cal column densities over a period of several years. Since 
the previous observations in 2002, the transverse velocity of 
the star has shifted the line-of-sight through the interstellar 
cloud by ~ 10 AU. 

From a consideration of the Doppler b parameter in each 
component, a rigorous upper limit was placed on the kinetic 
temperature of each component. Combining limits from each 
species, the temperatures of components A and B must be 
Tk,A < 67lt\v and Tk,B < 114^^4 K. However, the temper- 
atures may be lower if the turbulent velocity makes a sig- 
nificant contribution to the 6 values. The Cal/Kl ratio was 
used to derive logarithmic depletions of D(Ca)A ~ —1.82 
and D(Ca)B = -2.13. 

The Cal/Call ratio was used to derive estimates of the 
electron density in the two components. By assuming that 
electrons originate from photoionisation of neutral atomic 
carbon and assuming a value of 60% for the carbon deple- 
tion, lower limits on the density of ua ^ 7 x 10'^ cm~'^ and 
riB ^ 2 X 10* cm~'^ were derived. These values are far higher 
than typical densities in the diffuse ISM. Comparison of the 
observations with the chemical model of Bell et al. (2005) 
confirms the requirement for a high number density, with 
n = 5 X 10* cm~^ in the best-fitting model. 

The first detections of diffuse interstellar bands towards 
K Vel were made at two epochs, but only a loose upper 
limit of < 40% could be placed on the variation of DIB ab- 
sorption between the two epochs. The DIBs detected are 
weak given the measured _E(a_v), and an unusually high 
A5780/A5797 ratio of > 22 was found, amongst the highest 
known. This suggests that the carrier of A5797 is more sen- 
sitive to UV radiation than local density. The behaviour of 
DIBs towards k Vel appears to be similar to that observed 
in Orion by Cami et al. (1997). The weakness of the DIBs 
(per unit i?(B-v)) may be related to the high density of the 
material towards k Vel. 

Continued monitoring of interstellar absorption lines to- 
wards K Vel would be of interest, in order to determine the 
size of the structure from the time taken for the column den- 
sity to peak and subsequently decline. It would be beneficial 
to observe the C l fine structure lines in the ultraviolet to al- 
low derivation of the local density and kinetic temperature. 
Although this approach leads to some degeneracy in these 
parameters, this could be lifted via the existing constraints 
on temperature and density, though (nearly) simultaneous 
optical observations would be an advantage. A direct mea- 
surement of the H2 fraction would also be beneficial. 
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